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ABSTRACT 

The Muon g — 2 collaboration has measured the anomalous magnetic g value of the 
positive muon to within a relative uncertainty of 0.7parts per million. The result, 
a^+(expt) = 11 659 204(7) (5) x 10~^°, is in good agreement with the preceding data 
on and a^- and has about twice smaller uncertainty. The measurement tests 
standard model theory, which at the level of the present experimental uncertainty 
involves quantum electrodynamics, quantum chromodynamics, and electroweak in- 
teraction in significant ways. The analysis of the data on the anomalous magnetic 
g value of the negative muon is well underway. 
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1 Introduction 



The anomalous g values, a = {g — 2)/2, of leptons arise from quantum mechani- 
cal effects. Their precise measurement has historically played an important role in 
the development of particle theory. The anomalous magnetic g value of the elec- 
tron, tte, has been measured to within about four parts per billion (ppb) j2], and 
is thus among the most accurately known quantities in physics. Its value is de- 
scribed in terms of standard model (SM) field interactions, with nearly all of the 
measured value contributed by QED processes involving virtual photons, electrons, 
and positrons. Heavier particles contribute to Qe only at the level of the present 
experimental uncertainty. 

The anomalous magnetic g value of the muon, a^, is more sensitive than 
Oe to processes involving particles more massive than the electron, typically by a 
factor {m^/nieY ~ 4 • 10^. A series of three experiments jS] at CERN measured 
to within 7 parts per million (ppm), which is predominantly statistical. The CERN 
generation of experiments thus tested electron-muon universality and established 
the existence of a hadronic contribution to with a relative size of ~ 59 ppm. 
Electroweak processes are expected to contribute at the level of 1.3 ppm, as are 
many speculative extensions of the SM. 

The muon g — 2 experiment at Brookhaven National Laboratory (BNL) is 
conceptually similar to the last CERN experiment, and has determined of the 
positive muon with an uncertainty of 0.7 ppm from a sample of about 4 • 10^ decay 
positrons collected in the year 2000. The analysis of a^- of the negative muon from 
a similarly sized sample of decay electrons collected in 2001 is well underway. 

2 Experiment 

The measurements at the Brookhaven Alternating Gradient Synchrotron used a sec- 
ondary pion beamline. For most of the data taking periods, longitudinally polarized 
muons of about 3.1 GeV from forward decays were momentum-selected and injected 
into a 14.2 m diameter storage ring magnet |3| through a field-free inflector [S] re- 
gion in the magnet yoke. A pulsed magnetic kicker located at approximately one 
quarter turn from the inflector region produced a lOmrad deflection which placed 
the muons onto stored orbits. Pulsed electrostatic quadrupoles [7j provided vertical 
focusing. The magnetic dipole field of about 1.45 T was measured with an NMR 
system |H] relative to the free proton NMR frequency uop over most of the 9 cm di- 
ameter circular storage aperture. Twenty-four electromagnetic calorimeters |U] read 
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Figure 1: Top view of the g — 2 apparatus. The beam of longitudinally polarized 
muons enters the superferric storage ring magnet through a superconducting inflec- 
tor magnet located at 9 o'clock and circulates clockwise after being placed onto 
stored orbit with three pulsed kickers modules in the 12 o'clock region. Twenty- 
four lead scintillating-fiber calorimeters on the inner, open side of the C-shaped 
ring magnet are used to measure muon decay positrons (electrons). The central 
platform supports the power supplies for the four electrostatic quadrupoles and the 
kicker modules. 

out by 400 MHz custom waveform digitizers (WFD) were used on the open, inner 
side of the C-shaped ring magnet to measure the decay positrons and electrons. 
Muon decay violates parity, which in the laboratory frame results in a modulation 
of the number of positrons (electrons). 



above an energy threshold E. Here, Nq is a normalization, 7r ~ 64 /is is the dilated 
muon lifetime, A ~ 0.4 is an asymmetry factor, is a phase, and Ua is the angular 
difference frequency of muon spin precession and momentum rotation. 

The muon anomalous magnetic g value is evaluated from the ratio of the 
measured frequencies, R = Ua/ujp, according to: 



in which A = fi^/ fip is the ratio the muon and proton magnetic moments. The 
value with smallest stated uncertainty, A = fi^/ fJ^p = 3.183 345 39(10) jTU], results 
from measurements of the microwave spectrum of ground state muonium ^T] and 



N{t) = No{E) exp 





theory [Tg. 
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3 Data Analysis 



The proton NMR frequency uip and the muon spin precession frequency uja were 
analyzed independently by several groups within the collaboration. The values of 
R = uja/'^p and were evaluated only after each of the frequency analyses had been 
finalized; at no earlier stage were the absolute values of both frequencies known to 
any of the collaborators. 

3.1 The magnetic field 

During the data collection period from January to March 2000, a field trolley with 
17 NMR probes was used 2-3 times per week, 22 times in total, to measure the 
field throughout the muon storage region. Figure |2ti' shows the field value measured 
in the storage ring with the center trolley probe versus the azimuthal angle. The 
field is seen to be uniform to within about ±50ppm of its average value over the 
full azimuthal range, in particular also in the region near 350° where the infiector 
magnet is located. Since the field averaged over azimuth is uniform to within 1.5 ppm 
over the storage aperture (Fig. |2l3)5 the field integral encountered by the (analyzed) 
muons is rather insensitive to the precise location and profile of the beam. 

The probes inside the field trolley were calibrated with respect to each 
other during the data collection period using dedicated measurements in which a 
single NMR probe was plunged into the storage vacuum. This so-called plunging 
probe, as well as a subset of the trolley probes, were calibrated in situ with respect 
to a standard probe [T^ . 

The 22 measurements with the field trolley were used to relate the readings 
of 370 NMR fixed probes in the outer top and bottom walls of the storage vacuum 
chamber to the field values in the beam region, so that the fixed NMR probe readings 
could be used to interpolate the field when the field trolley was 'parked' in the storage 
vacuum just outside the beam region and muons circulated in the storage ring. 

For the data collection between January and March 2000, the field fre- 
quency Up weighted by the muon distribution was found to be, 

cjp/(27r) = 61 791 595(15) Hz (0.2 ppm), (3) 

where the uncertainty has a leading contribution from the calibration of the trolley 
probes and is thus predominantly systematic. The result was confirmed by a sec- 
ond, largely independent analysis, which made use of additional calibration data, a 
different selection of fixed NMR probes, and a different method to relate the trolley 



4 




Multipoles [ppm] 
normal skew 
Quad 0.24 0.29 
Sext -0.53 -1.06 
Octu -0.10 -0.15 
Decu 0.82 0.54 



azimuth [degree] radial distance [cm] 



Figure 2: The NMR frequency measured with the center trolley probe relative to 
a 61.74 MHz reference versus the azimuthal position in the storage ring (left), and 
(right) a 2-dimensional multipole expansion of the azimuthal average of the field 
measured with 15 trolley probes with respect to the central field value of 1.451 275 T. 
The multipole amplitudes are given at the storage ring aperture, which has a 4.5 cm 
radius as indicated by the circle. 

and fixed probe readings. The result from the data collection on the negative muon 
in the year 2001 is expected to have further improved uncertainty. 

3.2 The muon spin precession frequency 

About 4 ■ 10^ reconstructed positrons with energies greater than 2 GeV and times 
between 50 fis and 600 /is following the beam injection were available for analysis 
from the data collection between January and March 2000. Figure shows their 
time spectrum after corrections for the bunched time structure of the beam and for 
overlapping calorimeter pulses, so called pile-up [iHj, had been applied. 

The main characteristics of the spectrum are muon decay and spin preces- 
sion (Eq. H)), however, additional effects need to be considered as illustrated by the 
Fourier spectrum in Fig. Eb- These effects include detector gain and time instability, 
muon losses, and oscillations of the beam as a whole, so-called coherent betatron 
oscillations (CBO). The latter are caused by the injection of the beam through the 
relatively narrow 18(w) x 57(h) mm^ aperture of the 1.7m long infiector channel into 
the 90 mm diameter aperture of the storage region. Their frequencies are determined 
by the focusing index of the storage ring, and have been observed directly with fiber 
harp monitors that were plunged into the beam region for this purpose. Since the 
calorimeter acceptances vary with the muon decay position in the storage ring and 
with the momentum of the decay positron, CBO cause modulation of the time and 
energy spectra of the observed positrons. 

Numerically most important to the determination of Ua from the data col- 
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Figure 3: The time spectrum for 4 ■ 10^ positrons with energies greater than 2 GeV 
collected from January to March 2000, after corrections for pile- up and for the 
bunched time structure of the injected beam (left) were made, and (right) the Fourier 
transform of the time spectrum, in which muon decay and spin precession (cf. Eq. 2) 
has been suppressed to emphasize other effects. 



lected between January and March 2000 were the CBO in the horizontal plane, whose 
frequency was numerically close to twice the frequency uJa- When the modulations 
of the asymmetry and phase with frequency o^cbo.h — 2 x were not accounted for 
in the function fitted to the data, artificial shifts of up to 4ppm in the frequency 
values UJa determined from individual calorimeter spectra were observed. In the 
joined calorimeter spectrum, such shifts are largely canceled because of the circular 
symmetry of the experiment design. 

Several approaches were pursued in the analysis of Ua- In one approach, 
the time spectra from individual positron calorimeters was fitted in narrow energy 
intervals using a fit function as in Eq. ^ extended by the aforementioned number, 
asymmetry, and phase modulations. Other approaches made use of the cancellation 
in the joined calorimeter spectra and either fitted for the residual of the leading 
effects, or accounted for their neglect in a contribution to the systematic uncertainty. 
The results were found to agree, on uja to within the expected 0.5 ppm statistical 
variation resulting from the slightly different selection and treatment of the data in 
the respective analyses. The combined result was found to be, 

cua/(27r) = 229 07411(14)(7)Hz (0.7ppm), (4) 

in which the first uncertainty is statistical and the second is systematic. The above 
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Figure 4: Recent measuremens of and standard model evaluations using the 
estimates in Ref. j2S] of the lowest order contribution from hadronic vacuum polar- 
ization. 

frequency includes a correction of +0.76(3) ppm for the net contribution to the muon 
spin precession and momentum rotation caused by vertical beam oscillations and, for 
muons with 7 ^ 29.3, by horizontal electric fields P^. The systematic uncertainty 
has a leading contribution of 0.2 ppm caused by CBO. In the year 2001, an event 
sample of comparable size on the negative muon was collected. The storage ring 
was operated with two different values of the focusing index, which is expected to 
reduce the leading systematic uncertainty. 

4 Results and Discussion 

The value of was evaluated after the analyses of Up and Ua had been finalized, 

= 11659 204(7)(5) x 10"^° (0.7ppm), (5) 

where the first uncertainty is statistical and the second systematic. The result 
agrees well with the preceding measurements jHl HSl HHI and drives the present world 
average, 

a^(exp) = 11659 203(8) x 10"^° (0.7ppm), (6) 

in which the uncertainty accounts for known correlations between the systematic 
uncertainties in the measurements. Figure |3] shows our recent measurements of a^+, 
together with two SM evaluations discussed below. 

In the SM, the value of receives contributions from QED, hadronic, 
and electroweak processes, a^(SM) = a^(QED) + a^(had) + a^(weak). The QED 
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and weak contributions can, unlike the hadronic contribution, be evaluated per- 
turbatively, a^(QED) = 11658 470.57(29) x 10-^° [HI and a^(weak) = 15.4(2) x 
IQ-io jjg' ^ fjj-^g hadronic contribution is, in lowest order, related by dispersion the- 
ory to the hadron production cross sections measured in e~^e~ collisions and, under 
additional assumptions, to hadronic r-decay. Clearly, the hadronic contribution has 
a long history of values as new data appeared and analyses were refined. 

At the time of the PIC-2003 conference, a recent and detailed evaluation 
was the one by Davier and co-workers, which - unlike preceding analyses - incor- 
porated the low-energy e~^e~ annihilation cross section into hadrons by the CMD-2 
collaboration [201 ? e~^e~ measurements with improved accuracy in the 2-5 GeV 
energy region from BES, preliminary results from the final ALEPH analysis |22j of 
hadronic r-decay at LEPl, and additional data PS] from CLEO. Significant discrep- 
ancies between the e~^e~ and r data were found. 

The CMD-2 collaboration has since released a reanalysis of their cross sec- 
tion measurements and Davier and co-workers have provided updated estimates 
for the contribution to a^(SM) from lowest order hadronic vacuum polarization, 
a^(had, 1) = 696(7) x 10"^° from e+e~ data and a^(had, 1) = 711(6) x 10"^° from 
r data [23]. The authors refrain from averaging the values, noting that significant 
discrepancies remain in the underlying data in the center-of-mass region between 
0.85 and 1.0 GeV. Radiative return measurements at the e"'"e~ factories may reach 
the required precision to shed light on the situation, as might lattice calculations. 
Higher order contributions to a^(had) include higher order hadronic vacuum polar- 
ization [2n] and hadronic light-by-light scattering |27j . 

Open questions thus concern the SM value of a^, in particular its hadronic 
contribution, and the experimental value of a^- at sub-ppm precision. Stay tuned! 
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